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The optical conductance of a multiple scattering medium is defined as the total transmitted 
light of a normalized diffuse incoming beam. This quantity, analogous to the electronic conduc- 
tance, exhibits universal conductance fiuctuations. We perform a detailed diagrammatic analysis 
of these fluctuations. With a Kadanoff-Baym technique all the leading diagrams are systematically 
' generated. A cancellation of the short distance divergencies occurs, that yields a well behaved 

theory. The analytical form of the fluctuations is calculated and the theory is applied to optical 
systems. Absorption and internal reflections, which are taken into account, reduce the fluctuations 
significantly. 



o 

CD 

Q 

00 



X 



I. INTRODUCTION 



> 

' When studying the transmission of optical waves through random media one may consider fluctuations in the 
\l optical conductance. Analogous to the electronic conductance, this quantity is obtained when, first, the incoming 
beam is monochromatic and diffuse, and second, all outgoing light is collected. Since the electronic conductance 
is known to exhibit so called universal conductance fluctuations (UCF), the same fluctuations are expected in the 
optical conductance. Measurements of this quantity would constitute a cornerstone in the analogy between optical 
\ and electronic mesoscopic systems. 
"Y"^ ■ We briefly review the situation for mesoscopic electron systems The electronic conductance of mesoscopic 

' samples is known to show reproducible sample to sample fluctuations. Since the fluctuations are a consequence of 
, scattering from static impurities, they are static. Their magnitude is independent of the sample parameters such as 
I ' the mean free path, the sample thickness and the average conductance. Hence they are called universal conductance 
' ^ ' fluctuations. The mean conductance in the considered regime comes from multiple scattered diffuse electrons. The 
UCF are a consequence of interference of multiple scattered waves, causing correlations between two diffuse paths. 
^ ' Therefore the fluctuations are much larger then one would obtain classically by modeling the system by a random 
resistor array, in which interference effects are neglected. That approach is valid only on a length scale exceeding the 
phase coherence length, where the fluctuations reduce to their classical value. The conductance fluctuates when the 
phases of the waves in the dominant paths are changed. This happens, of course, if one changes the position of the 
$^ ' scatterers, eq. by taking another sample. One may also keep the scatterers fixed but apply a magnetic field or vary 
the Fermi energy. (In optical system one can vary the frequency of the light.) In all these cases one modifies the 
phases of scattered waves, so that different propagation paths become dominant. 

Transport through mesoscopic systems is not only studied in electron systems, but also using multiple scattered 
classical waves such as sound, microwaves, and, particularly, light. The origin of the mesoscopic phenomena is the 
interference of multiple scattered waves and in first approximation all systems are described by the same equations, 
namely the scalar wave equation, see for instance van Haeringen et al. [|6| for a review. Thus one expects that similar 
large fiuctuations in optical systems are present in the diffuse transmission regime. An advantage of optical systems 
over electronic ones is that optical systems are much cleaner: no equivalents of phonons or electron-electron interactions 
are present. Indeed very accurate measurements of the enhanced backscatter cone 0, correlation functions |||,^ and 
intensity distributions were performed. Nevertheless, to the best of our knowledge, the optical analog of the UCF 
has not yet been observed in optical systems. Such experiments turn out to be difficult. Although the magnitude of 
the fiuctuations is universal, they occur on a background of order where g is the dimensionless conductance (in 
optical experiments one typically has g ^ 10^). The relative value of the fluctuations to the background is thus 1/g, 
so that the C3 correlation function is of order 1/g^, typically of order 10~^. For electrons this problem is absent 
as moderate values of g are achievable. This is also the reason that electrons are easier brought near Anderson 
localization, for which g has to take a critical value of order unity. In the electronic case the moderate values of g, 
combined with very sensitive techniques for current measurements, have led to many observations of the universal 
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conductance fluctuations. Recent optical experiments suggest, however, that the optical analog of UCF should just 
be experimentally accessible, as the measurement of the third cumulant in the total transmission was reported [[lO| . 
This quantity is of the same order as the optical UCF, namely It is expected that similar techniques can be 

applied to measure the optical UCF. Microwave scattering is also interesting as it combines lower values of g with 
many of the advantages of optical systems. 

We are not only interested in just the size of the fluctuations, but also in the somewhat more general frequency 
correlation. A change in the frequency alters the interference pattern, just as occurs by changing the magnetic 
field or the Fermi energy in the electronic case. It is known from experiments that using the frequency as tunable 
parameter provides a good way for measuring the fluctuations. In contrast to electronic systems, in optical systems 
both angular resolved and angular integrated measurements of the transmission are possible. Therefore various 
transmission quantities can be measured in optical systems, each with its particular frequency correlation function. 
In a recent review by Berkovits and Feng the different correlations and their physical interpretation are discussed 
extensively. Denoting the transmission from incoming channel o (wave coming in under angles 9a, ( j^a) to outgoing 
channel b (waves transmitted into angles 9b, cj)^) as Tab, the correlation functions can be classifled as 



J.,(c.)r ,(c. + Aa;)) ^ ^ ^ cr\^u.) + cr^iAco) + Cf^'{Au;). (1) 

{J-ab[^)){-l-cd[^ + ^^)) 



The unity just comes from the product of averages. The Ci-term in the correlation function is the most important 
one if the angular resolved transmission Tab is measured. The reason is that Ci is of order unity when both the 
incoming directions a and c, and the outgoing directions b and d are, pairwise, close to each other. When dealing with 
one single monochromatic plane wave (a = c), Ci describes the correlation of the bright and dark speckle pattern. 
Diagrammatically, Ci is the sum of all reducible diagrams, that is to say, just as the unit contribution in Eq. (Q), it 
is equal to the product of two averages. 

If instead of measuring light in one outgoing channel, all outgoing light is integrated, the sharply peaked and short 
ranged Ci correlation function is overwhelmed. By collecting the outgoing light, the total transmission To = X^b^afc 
is measured; experimentally this is commonly done using an integrating sphere In this setup the C2 correlation 
function, which has a much smaller peak value but is long ranged, contributes for all outgoing angles and becomes 
dominant. Its long range character arises because, due to interference of the diffuse light paths, the outgoing amplitudes 
are pairwise in phase. The C2 correlation, which still depends on the angles of the incoming beams a and c, is of 
order g~^. The C2 corresponds to a diagram where the two incoming diffusons interact through a Hikami- vertex. 

Finally, the C3 term is dominant when the incoming beam is diffuse, and all outgoing light is collected, so that, 
just as in electronic systems, the conductance g ^ T = '^ab'^'^b is measured. In that measurement contributions 
where a and c are far apart are dominant. In contrast to previous case, now also these incoming amplitudes must be 
pairwise in phase. This occurs in a diagram where the two incoming diffusons interact twice, so that a loop occurs; 
in principle, further loop insertions to it also contribute. In spite of the fact that C3 is of order g^^, it is dominates 
over the Ci and C2 terms as it has contributions for all incoming and outgoing angles. 



The Ci and C2 correlation have been studied in detail, both experimentally [y,y_3 14 1 and theoretically [|5|. It was 
also shown in experiments [ p^ and in theory that absorption and internal reflection, neglected in the earliest 

calculations, significantly reduce the correlations. Among other methods, the Ci and C2 correlation were successfully 
[ p5| , p^ calculated using a diagrammatic technique based on the Landauer approach (l^j2^ . One might hope that the 
calculation of the C3 or UCF in this approach is also straightforward. It is well known, however, that the calculation 
in the Landauer approach is quite cumbersome, since divergencies show up on scales of one mean free path when the 
problem is treated on a macroscopic level using diffusons. 

In order to circumvent these difficulties, one is tempted to use the Kubo approach, often used in electronical systems 
to calculate the UCFJ&^. Furthermore, the results for the conductance obtained by Kubo or Landauer formalism 
should be identical pl| , |22| . Yet the Kubo approach can not be applied directly to optical systems, since it is not clear 
how external lines should replace current vertices, and how absorption and internal refiections are to be included. 
Therefore, we use the Landauer approach. 

Technically, the difficulties in the Landauer approach are caused by the vertices for partner exchange of two diffusons, 
the so-called Hikami boxes. Each Hikami box brings the square of the internal momentum, whereas the current vertices 
are momentum- independent in the Kubo formula. As a result, the integral over the internal momentum of the closed 
loop is convergent in the Kubo formula, while naively divergent in Landauer approach. 



Two studies of the C3 in the Landauer approach are known to us. In the first, Kane, Serota and Lee |23| consider 
electronic systems and make elegantly use of current conservation to derive an expression for the correlation function. 
Although in optical systems the conserved quantity is not the intensity but the energy, their prediction applies to 
optical systems as well, since it amounts to a result for the same sums of scattering diagrams, involving different 
parameters only. This result has not been confirmed by a direct derivation, however. Moreover, since it relies on 
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a conservation law, it is not clear what happens when absorption is present. The second study was performed by 
Berkovits and Feng After giving a very clear discussion of the problem, these authors calculated one of the 

macroscopic diagrams (presented earlier by Feng, Kane, Lee and Stone and subtracted the divergent parts by 
hand. In this way the correct order of magnitude and the qualitative frequency dependence was obtained. 

It is our goal to clarify the situation by calculating the optical C3 diagrammatically. A complete analysis of all 
leading diagrams is needed, in which finally the divergent parts should cancel. We specialize to the case where the 
conductance is measured. In this setup the amplitudes of the incoming and outgoing diffusons are exactly in phase. If 
the C3 correlation is measured as a (small) part of the correlation in the angular transmission or total transmission, 
this phase condition need not be full filled. Other contributions of order Xjg^ which are angular dependent are then 
present. These contributions are both the diagrams presented below, with different decay rates for the incoming 
and outgoing diffusons, but also new diagrams contribute. In conductance measurements these complications do not 
occur. Our fundamental approach immediately allows for inclusion of effects due to boundary layers and absorption. 
Our calculations, although specialized to optical systems, are valid for any mesoscopic system. 

Tacitly we assumed that higher order correlations can be neglected. Altshuler et al. ||2^ have shown that for 
small n the rtth cumulants scale as so that for g ^ 1 the fluctuations have indeed a Gaussian distribution; 

the far tail of the distribution is predicted to be log-normal. It would be interesting, however, to determine the full 
distribution function of the fluctuations, as was done in experiment for the angular resolved transmission |l3| and the 
total transmission [lO| ] , and recently in theoretical work by Nieuwenhuizen and Van Rossum . 

The outline of this paper is as follows. First, we introduce the basics of the diagrammatic technique and describe 
diffuse transport of light. In section [II we present the long distance diagrams and analyze the divergencies arising 
from these diagrams in the diffusion approximation in IV. Next, in section ^ we develop a Kadanoff-Baym theory 
in order to generate all relevant scattering diagrams. In section ^ we show that the diver genc ies indeed cancel if all 
diagrams are analyzed in detail. The general form of the fluctuations is calculated in section VII and applied to optical 
systems, where absorption and internal refiections may be present. The presentation is closed with a discussion. 



II. DIFFUSE TRANSPORT OF LIGHT 



In this section we discuss standard aspects of diffuse light transport, such as the amplitude Greens function, the 
diffuse intensity and the transport equation. We shall consider the situation of point scatterers in a medium that has a 
dielectric constant different from its surroundings. We employ the notation and results of the paper by Nieuwenhuizen 
and Luck The main results of the present section are the expression ( ^7|) for the C3 correlation function, and 
the expression ( ^9| ) for the universal conductance fluctuations. These expressions involve diagrams with incoming 
total-flux diffusons Cin defined in (^8|), outgoing total-flux diffusons £out defined in (^6|), and internal diffusons Cmt 
defined in ( p5| ) or (|3^) . Readers not interested in the microscopic background of these relations may skip the details 
of the present section. 

We consider a quasi-one, quasi-two, or three dimensional slab of thickness L and area A = W'^'\ with {W > L). 
which contains static, isotropic point scatterers. As usual a scalar approximation is made for the electromagnetic field 
of the light For bulk properties this is justified since the polarization is scrambled after a few scattering events. 
For other applications, such as acoustic waves and spinless electrons, the scalar property is immediate. The scalar 
wave equation at given frequency uj reads 

V2V^(r) + ^e(r)V(r) = 0, (2) 

where c is the vacuum speed of light. We shall consider a medium with dielectric constant eo and density n of small 
spheres with dielectric constant 62 and radius oq, located at random positions Ri. Going to the limit of point scatterers 
we get 

e(r) = eo + ^7rao(e2 - £0) J^^^^^^i^l < z < L, (3a) 

i 

= €1 z < 0; z > L, (3b) 

where ei is the dielectric constant in the surrounding medium. The wave numbers in the surrounding medium and in 
the random medium are. 
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respectively. In the bulk the average retarded Greens function (amplitude Greens function) in three dimensions with 
momentum p reads 



Gip)^{p^-e-nt) \ (5) 
For small scatterers the scattering becomes isotropic with an effective strength described by the i-matrix |Q 

^TZ ' (6) 



u uq Att 

where u = 47rfc^ao(e2 — eo)/3 is the bare scattering strength and uq is an internal parameter of the point scatterer. 
The mean free path is defined as 

47r 

i--^, (7) 
nil 

in which n is the scatterer density. The metallic mesoscopic regime is defined by: k£ ^ 1 and £ L . Eq.(||) satisfies 
the optical theorem Im(i) — kit/ Ana, where a is the albedo of the scatterer. The case of no absorption, which we 
treat first, corresponds to a = 1. In electronic systems absorption does not take place, but a similar behavior arises 
from de-phasing; in that case the phase coherence length plays the role of absorption length and, for our purpose, has 
to be larger than the system size. 

For disordered electron systems one deals with weak s— wave scattering and the diagrammatic expansion can be 
carried out in second order Born approximation. In experiments on optical systems however, efhcient scattering is 
achieved by taking resonant scatterers and strictly one now has to calculate the full Born series |^8| . In the main part 
of this work we study how the cancellation of some short range contributions takes place. This we do in the second 
order Born approximation, where t « u + iu^k / At: and £ ~ Air/nu'^. Working with the full Born series would introduce 
a large number of extra classes of diagrams so that the total number of diagrams increases dramatically. We expect, 
however, that the cancellation of divergencies, when shown in the second order Born approximation, also holds if the 
full Born series is considered. It is widely expected and confirmed in one of our earlier papers that for static quantities 
the difference between the two should only be a renormalization of the mean free path in the final results |2^. We 
thus continue to work with the full Born series where possible. The cancellation of short range divergencies will only 
be shown within the second order Born approximation. 

Intensity transport in transmission is dominated by ladder diagrams or diffusons. A diffuson is made up by pairing 
one retarded and one advanced propagator sharing the same path through the sample, as is depicted in Fig. |^. 

Although the transport deep in the bulk is accurately described with the diffusion equation, the precise behavior 
near the (reflecting) surface has to be derived from the Schwarzschild-Milne integral-equation Consider a plane 
wave with unit flux impinging on the sample under angle 9a 

y'Aki COS Oa 

where p = {x,y) is the transversal coordinate and Qa = fci sin0a(cos(/)a, sini/ia) is the two dimensional transverse 
momentum of the incoming beam. Inside the slab the unscattered part of the intensity decays over one mean free 
path. This is described by (||) since its pole lies slightly away from the real axis, due to the imaginary part of t. The 
source of diffuse intensity is the first scattered intensity. 

Saiz) = T{pa) ntt |^i"(r)p = MMe"^/^'^", (9) 

Akiflat 

where /ia = cos 6'^ involves the angle 9'^ of the refracted beam with respect to the z-axis. The source is proportional 
to the intensity transmission coefficient of the boundary between the two dielectrics 



, , AuJ 11^ - 1 + llm? , , 



j^212 



The twice scatterered intensity follows from Sa as ntt J dr'\G{r,r')\'^ Sa{r'). The geometric sum of first, second, 
third,..., times scattered intensity is the multiple scattered or diffuse intensity Ca- It is generated from the first 
scattered intensity by the Schwarzschild-Milne equation 
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Ca{z) = Sa{z) + r dz'M{z, z')L^{z'). (11) 

Jo 

The kernel M reads 

z') = Mb{z - z') + Ml{z + z') + Mi(2L - z - z'), (12) 

which contains the bulk term 

Mb(z-z')= r J^e-l^-^'l/''^ (13) 
and layer terms describing internal reflections at the interfaces at z = and z = L. The layer term reads 

'^^M^i rr^l ..\\^-(z+z')/i,t 



Ml{z + z')= f ^(l_r(^))e- 



(14) 



It involves the total path length of the radiation going from depth z to the boundary at z = 0, reflected there and 
going to depth z' . In the bulk of the slab (several mean free paths away from the boundaries) the diffuse intensity 
has a simple behavior 

^inj-^^ ^ 47rr(^a)Ti(pta) L + Zp - Z ^^^^ 

° ki£m Afia L + 2zo 

Here ti is a factor that determines the limit intensity of a semi-infinite slab; of course, it depends on the incident angle. 
The length zq is called the "injection depth". For isotropic scattering one has [|2^ zq — 0.7104£. However, it becomes 
larg er when there is a mismatch of the indices of refraction between the scattering medium and the surroundings 
|p8| . Note that (^5|) satisfies the diffusion equation V^/CJf = 0; the more complicated Schwarzschild-Milne equation is 
needed to fix the parameters of its solution. 

In the electronic conductance measurements, however, waves coming from all directions are involved. In the 
Landauer formula for the conductance, 

G^^Y^^ab, (16) 
ab 

one needs to sum the transmission coefficients Tab of waves with unit flux coming in channel a and going out to 
channel b. In an optical experiment, however, the incoming diffuse beam may have angular weights that differ by a 
factor of order unity. When using integrating spheres one measures the outgoing intensity, rather than the outgoing 
flux. As compared to the electronic case, it does not bring a factor fib to the weight of the outgoing channel. As these 
differences between electronic and optical measurements only lead to different numerical pre-factors, we calculate the 
optical UCF with the same weights as in the electronic case. Summing (|^) over the channels a yields a source for the 
diffuse intensity 

S{z) = Y,Sa = ^ TrfMe-/^^". (17) 
* Jo 



This is again the input in the Schwarzschild-Milne equation(|ll|). In the bulk the intensity now has the diffusive 
behavior 

This object has been termed the incoming total- flux diffuson [p5| . It has the same depth dependence as C°^^, but 
contains a different prefactor. Deflning 



" kfm^Afia 

we can verify the sum rule 



(19) 
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Here we used that k = mki, the definition /i = -^/l — Q'^/k^ and (2.30) of Ref. [2^. One thus has 

£r„(z) = eaA„(^). (21) 

On the outgoing side, radiation emitted at a point r = {p, z) inside the slab will propagate to a point {Z, p') outside 
the sample {Z > L) as described by the Greens function of semi-infinite medium with dielectric function e(r) = eo for 
z < L and e(r) — ei iov z > L 

G{p, z;p',Z) = jYl G{z- Z- Q)e^Q^P-p'\ (22) 
Q 

in which G{z; Z; Q) is the 1-d Fourier transform of Eq. (|^), 



G{Z-Z-Q)^- ^P{L-z)+^p{Z-L) ^ P=yp~QM^, p=./fc2_Q2_ (23) 

In the far field {Z » L) the total transmitted intensity reads 

/ dV'|G(p, z; p', Z)\' = jY. (24) 

Q 

Since = fci cos^q = kpa^ the according flux is 



(25) 



It leads to a source S{z) = 47r4>(z)/£ in the Schwarzschild-Milne equation. The outgoing total-flux diffuson in 
transmission therefore reads 

r / \ k z + zo 



Apart from a reflection, this expression differs by a factor 4 from (18). The pre- factors would be the same if our 
Greens functions were multiplied by a factor 2; this amounts to the same as taking a kinetic term V2/2 rather than 
V^, such as occurs in electronics in units where Ti = m = 1. 

For outgoing intensity in direction h one has, in analogy with (|2l|). 



£^,t(z) = eb£out(^) (27) 

We call the incoming and outgoing diffusons, external diffusons, because they are connected to the outside of sample. 
This in contrast to the internal diffusons that begin and end at interference vertices inside the medium. Away from 
the surface the internal diffusons obey the well known diffusion equation: V2£mf(r) = 127r5(r — r')/i^ . 
The result for the angle-resolved transmission of Nieuwenhuizen and Luck pGf can be written as 

{T)ab = eaeb{T), (28) 

where (T) is the average conductance 

in dimensionless units. Restoring units for the electronic case one has for the average conductance 

All the above can be generalized to include absorbing scatterers, frequency differences between incoming beams. For 
an electronic system the correspondent effects would be a change in Fermi energy, rather than a change in frequency. 
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and the effect of a finite incoherence length, rather than the effect of absorption. For this more general case the 
diffusion equation reads 

[-V^ + k' + in]C,nt{r) = ^5{r - r'). (31) 

The inverse absorption length k is related to the albedo a, as = 3(1 — a)/£^; — Aoj /D is the ratio of the frequency 
difference of the amplitude propagators and the diffusion constant D. It holds that D = ^vsi, where ve is transport 
speed poj ]. In the slab geometry after Fourier-transform in the z-direction the diffusons obey: 

- ^ Ant(^, z; M) -f M^C,nt{z, z'- M) = ^5{z - z'), (32) 
where we have defined the decay rate 

= Q'^ + K,'^ + in. (33) 

The complex parameter M describes the exponential decay C ~ exp(— Mz) of the diffuse intensity in the z-direction. 
The solutions to the diffusion equation are a linear combination of hyperbolic sines and cosines. In conductance 
measurements using integrating spheres external diffusons with momentum or frequency terms yield no contribution, 
and thus M — k for the external diffusons, i.e. 



4fc sinh k(L — z) -|- kzq cosh k.{L 
I (1 -|- k^Zq) sinh kL -|- 2kzo cosh kL ' 



^ fc sinh Kz -I- KZo cosh kz (34b) 

°" (1 + K^^^) sinh kL -|- 2kzo cosh kL 

Whereas for the internal diffusons 

, 127r [sinh Mz + Mzo cosh Mz] [sinh Af(i-z')-hMzo cosh M(L-z')] 

A„t(z, z ; M) = — (A^ + M^zl) sinh ML + 2M^zo cosh ML ' ^^^^ 

where it is assumed that z < z', otherwise z and z' must be interchanged on the r.h.s. With equal indices of refraction 
inside and outside the sample the extrapolation length zq is, as stated above, 0.7104i! and thus the terms involving zg 
yield contributions of the order £/L. For optically thick samples {L ^ £) one has 

^ , , 127r sinh Mz sinh M (L - z') 

^-(^'^ = ^ m^uLaTl ^''^ 

Consistent with this expression the boundary conditions become approximately £i„t(0, z') = Cint{L, z') = 0. If there 
is an index mismatch, the surfaces partially reflect and the extrapolation length zq increases [|2^ . Internal reflection 
becomes especially important when zq becomes comparable to the sample thickness. Such may occur for large index 
mismatch and moderate thicknesses. 

The C3 correlation function, defined in (^, involves incoming diffusons ^ and outgoing ones >Couf Due to the 
factorization of external direction dependence, see (|2|), (p7|), and (|2^), it cancels from C3. We can write 



where Q is the two-dimensional transversal momentum. The function F is the main object to be determined in this 
paper. It is thus calculated at fixed Q and with external diffusons being total-flux diffusons. One finds from (|l]) and 
( p7| ) the conductance fluctuations 

CTiK,n) = {T {u;)T {u; + Aco)) - (T(w))(T(u; -t- Aw)) = ^ (T)afa(T),dC3 

abed 

= J2FiQ,n,n) 

Q 

= F(0, K, n) quasi 1-d, (38a) 

= W j '^FiQ, ^) quasi 2-d, (38b) 

= W^ j ^F{Q,n,n) 3-d. (38c) 
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For electronic systems one finds 

(G(fc)G(fc + yn)) - (G(fc))(G(fc + ^m)) = (^^^ CT{o,n). (39) 

These results can be extended for other geometries. If the width is comparable to the thickness of the slab, the 
momentum integral discretizes into a sum over transversal eigenmodes The result can be generalized further to 
arbitrary geometries by taking x, y dependence into account, and calculating the diffusons using appropriate boundary 
conditions. 



III. LONG RANGE DIAGRAMS 

At this point we address the structure of the leading diagrams for the correlation functions, defined in (|^). For 
all diagrams there are two incoming advanced fields, which we momentarily term i and j, and two retarded ones, i* 
and i* . The first term on the r.h.s. in Eq. (|^) follows from the diagram where i and i* are paired into an incoming 
difFuson, and the same for j and j* . These diffusons have no common scatterers, so for this contribution the expression 
factorizcs into a product of averages. For the Gi term in ([T]) such a factorization also takes place jl2j. However, in 
this term the pairings are ij* and ji* . In the G2 correlation function there arc two terms. In the first the incoming 
diffusons have pairings ij* and ji* . These diffusons interfere in some point in space, where they exchange partners. 
The outgoing pairings are then ii* ^ jj*. It is this term that contributes in measurements of the total transmission. 
The time-inversed diagram, with in- and outgoing part of the diagram interchanged, also contributes to Eq. (|^). 
However, it does not contribute when this expression is summed over b and d, i.e. in total transmission. 

Interference of diffusons occurs if they exchange an amplitude in the presence of a common scatterer. In a dia- 
grammatic language this is described by so-called Hikami- vertices [^,^. They are depicted as shaded polygons, the 
four point vertex represents the interference of four diffusons, the six point vertex connects six diffusons. The shaded 
polygons in the figure indicate that the vertices are dressed, in second order Born approximation the dressed four 
point vertex is the sum of three diagrams, the six point vertex is the sum of sixteen diagrams |]3^ . 

In electronic experiments where one measures the conductivity, and in optical experiments where one uses an inte- 
grating sphere on the incoming side for creating a diffuse beam, the two amplitudes in the diffuson must have exactly 
the same phase. Therefore, the incoming diffusons cannot have a momentum or frequency difference, and the pairing 
must be ii* and jj* . In Fig. |^(a) the incoming diffusons interfere somewhere in the slab. In a diagrammatic language 
the diffusons interchange a propagator so that the pairing is changed into ij* and ji*. Propagation continues with 
these diffusons, which, due to the different pairing can have non-zero frequency difference and non-zero momentum. 
But in order to be dominant, the outgoing diffusons can also not have a momentum or frequency difference. Therefore, 
somewhere else in the slab a second interference occurs. Again exchanging an amplitude, the original pairing, ii* and 
jj*, is restored and the two diffusons propagate out, see Fig. H(a)i. Some other contributions occur as well. Whereas 
the incoming and outgoing pairings are always ii* and jj*, the internal ones may be different. In Fig. |^(a)ii the first 
incoming diffuson meets an outgoing diffuson and amplitudes are exchanged. These internal diffusion lines meet at 
a second point where the original pairings are restored. It is clear that in this process the intermediate paths are 
traversed in time reversed order. Due to time-reversal symmetry they give a similar contribution as previous Fig. H(a)i. 
In Fig. ^(b)i, for instance, a diffuson breaks up such that one of its amplitudes makes a large detour, returns to the 
breaking point and recombines into an outgoing diffuson. The second incoming diffuson crosses this excess path of 
the amplitude, and one of its amplitudes follows exactly the same contour as the one of the first diffuson. The fourth 
amplitude resides and finally recombines with its original partner amplitude to form an outgoing diffuson. Finally, 
Fig. ^(c) depicts the situation where only one internal diffuson occurs. Its endpoints must lie within a distance of 
a few mean free paths. Because of its local character this class does not show up in the final result; it is needed, 
however, in the regularization process, since it contains terms that cancel divergencies from the other two classes. 

Due to the diffusive behavior of the internal propagators we call all these the "long range diagrams" . Their internal 
lines are diffusons (ladder diagrams), that is to say, in these lines there can be an arbitrary number of scatterers. Note 
that these long range diagrams also include terms with only a small number of scatterers, e.g. one or two. The latter 
contributions are of course not really long range; however they contribute to the geometric series that represents the 
ladder diagram. We state this explicitly, since below we will discuss some unexpected problems of these short ranged 
contributions to the long range diagrams. There are also some special short ranged contributions. Due to various 



subtleties, they resist a general treatment; their calculation is postponed to section VI. 

Another class of diagrams can be constructed by taking the upper half of Fig. P(a)ii and combining it with two 
lower half of Fig. |^(b)i. This diagram contains again two four-point vertices and also two internal diffusons. It is 
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easily seen that incoming amplitudes ii* and jj* are mixed into ij* and ji*. These new diagrams are 1 / ^^-corrections 
to the angular dependent C2 and as such do not contribute in conductance measurements. 



IV. DIVERGENCIES IN THE DIFFUSION APPROXIMATION 



The diagrams for the conductance fluctuations contain a loop; the two internal difFusons have a free momentum, over 
which one has to integrate. In Fig. ^(a) this momentum is denoted q. Physically, one expects important contributions 
to the conductance fluctuations if the distance between the two interferences vertices ranges from the mean free path 
to the sample size. But in this section we will show that the q— integral for the long range diagrams diverges for 
large momentum, i.e. when the two interference processes are close to each other. The standard picture of diffuse 
transport with diffusons and interference described by Hikami vertices, that works so well for loop-less diagrams such 
as the C2 correlation function ||l5[], and the third cumulant of the total transmission now becomes spoiled by 



these divergencies. In section VI wc solve the problem by going back to mesoscopic scales, and consider all scattering 
events. 

The problem becomes clear if we calculate the diagrams of Fig. First, we need the expressions for the Hikami 
vertices. In order to derive the vertices, the momenta are expanded to leading order in (q£). For the large qi this is 
in principle not allowed, but in practice it could still work. The four point vertex and the six point vertex are found 
by summation of the bare vertex and its dressings. The calculation has been reported several times in the literature 



Ll|, P9y32| . We include the effects of absorption and frequency differences and find 

4 

2 



Hi = /i4[-qi • qs - q2 • q4 + ^ J^^cn^ + + (40a) 

He = -K[c{i •q2 + q2 -qs + Qs •q4 + q4 -qs + qs -qe + qe -qi + ^{^i + -^i^l + T^^i)]- (40b) 



2*2 

We call these the "invariant" forms of Hi and iJg as they are unchanged under the shifts q^ ^ q^ — \ Yl']=i^j^ 
qi — ^ qi — g qj' respectively. We defined the pre-factors as 

= 48^' = 96^- 

The momenta of the diffusons that are attached to these vertices are denoted by q^ , where the diffusons are numbered 
clockwise on the vertex and their momenta are directed towards the vertex. In the actual calculations the Fourier 
transforms in the z-direction of the vertices are used. Compared to previous results of, for instance Hikami p^ , the 
vertices contain additional frequency and absorption terms. According to the diffusion equation (^) these extra terms 
together with the q^ terms, lead to a source 5{z — z'). For external diffusons, such terms are neglected as they bring 
contributions of the order £/L. This approximation simplified the calculation of for instance the long range correlation 
function [p^ . For the internal diffusons, however, the source terms are of leading order and cause divergencies. They 
correspond to the situation where the two interferences take place within a distance of a few mean free paths. 

As an example we calculate the diagram presented in Fig. ||(a)i. This diagram was first depicted by Feng, Kane, 
Lee and Stone and considered in detail by Berkovits and Feng ||ll[]. These authors pointed out that a short 
distance divergency appears. For the case of external momenta approximately zero, the Hikami-box ( ^0|) yields 
i?4(q, 0, — q, 0) = 2hiq^, while the internal diffuson has the form Cint{q) = 127r/(£"^(7^). Omitting the external lines 
the diagram Fig. H(a)i. then simply leads to 

4 



77|(q,0,-q,0)/:L(9) = 44 / 7^9° = 4<5''' (r = 0), (42) 



(2^)3 "4 vm:", ^' ^ 4fc4 J (27r)3^ 4fc4 

which is indeed a cubic divergency in three dimensions, and, more generally, a c?-dimensional divergency in d dimen- 
sions. As it is arising from the physically innocent situation where the two interference vertices are close to each 
other, we expect that the divergency has to disappear finally. 

We now calculate the diagram for the slab geometry. For simplicity we first consider a quasi one-dimensional system 
in which frequency differences and absorption are absent, therefore the decay rate ("mass") vanishes, i.e. M — for 
all diffusons. (Beyond quasi Id one would have to take non-zero M = Q and sum over the allowed Q.) From Fig. ||(a)i 
one directly reads off its corresponding expression Fa.i 
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F ■ — 

^ a. I 



(43) 



We label the two incoming diffusons 1 and 3, the outgoing ones 2 and 4, and the internal ones £(2:5, zy) and C{zq, zg), 
(with 5,7 at z and 6, 8 at z'). The real space expressions for the Hikami boxes become 



H^iz') = ft4[5.,5., + lid., + d.,f - dl - dU, 



(44a) 
(44b) 



in which dz^ is the derivative of the corresponding diffuson; after performing the differentiation ^1,3,5^ should be put 
equal to z, while 22,4,6,8 should be put equal z' . Keeping 22,4,6,8 fixed, we obtain for the z-integral after some partial 
integrations 



dzHi{z)C\n{zi)Ci-n{z^)Cint{zZi Zq) Cint{z-; , Zg) 



dz[Cint{z, Z6)6{z- Zs) + C^nt{z, Zfi)S{z- Z(i)\C-^^{z) + 2hi / dzCi„t{z, Ze)Cint{z, zs)C[^{z) 



Cintizs, ZG)Cf^{zs) + 



127rft,4 



A«t(26,28)An(^6) + 2/l4 / dzCtntiz,ZG)Cintiz,Z8)Ci^{z). 



Here we also used the diffuson equation, which in this simplified case reads d^Cin = and d'^£int{z, z') — \2ti5{2 
z')/l^. Also carrying out the z'-integral we find after performing again some partial integrations 



dz' Hi{z')Cont{z2)Cout{zA) 



127r/i4 



Ant(28, Z6)A^n(^8) + 



£3 



Cint{z6,Zs)Ci^{ze) +2h4 / dz£int{z, Z6)£int{z, zs)£[n{z) 



4fc 

I 

+- 



J{0) I dzCl{z)Cl,,{z) 



dzC„u{z,z) +£2^(^)£;2^^(^) + 



+4/1I j dz' j dzCi^^{z,z')C'^{z)C'^^,{z'). 



The spatial derivative of C{z) is denoted C . All diffusons are simple linear functions in this case, yielding 

2 14 



(45) 



(46) 



Note that the prefactors of the diffusons and the Hikami boxes have canceled precisely. This is closely related to the 
universal character of conductance fluctuations in electronic systems, see (39). 

The term (5(0) is a linear divergency, which is the cause of all troubles. In the three-dimensional case one has to 
take Q ^ Q. The 5{Q)L term will occur also for transversal momentum Q ^ 0, so that the Q-sum yields the cubic 
divergency 



(2^)2 



(47) 



as expected from the above bulk consideration. 

We now give the results of all diagrams of Fig. ||. We no longer restrict to the M = case. The expressions are 
labeled according to the diagrams in figure. Fa, Fj, and Fc. In the diagrams of Fig. ^(a), the decay rates of Eq.(|3^) 
for the internal diffusons are each other complex conjugate, M and M*. In the diagrams of Fig. g(b) both internal 
diffusons have the same decay rate. Using the definition of the Hikami vertices and the diffusion equation, we obtain 
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Fa{M) 



Fb{M) 



Fc{M) 




(48b) 



(48a) 



(48c) 



where we used the short hand notation that in the single integrals all incoming and outgoing diffusons have argument 
z. To obtain the variance of the fluctuations in quasi-one dimension, _F's is evaluated at transverse momentum 
(5 = 0. The internal momentum Q enters the equations via the decay rate of the internal diffusons, defined by 
NP = Q-^ + K-^ + zil. In two and three dimensions a sum or integral over the transversal momentum has to be 
performed. 

In (^8|a,b) we can distinguish three contributions to F . In the first term, both Hikami boxes operate on the internal 
diffusons, yielding in the diffusion approximation a delta function evaluated in zero. The resulting term is independent 
of all the momenta of the external diffusons. It can be seen from the diffusion equation that a diffuson decays rapidly 
if its momentum becomes large. Terms of the diffusons with few scatterers are dominant at large momentum; they 
cause the divergence. Our present description of these processes is incomplete. In order to see the cancellation of this 
divergency, calculation of the long range diagrams is not sufficient, so that the short distance processes have to be 
examined in detail. This will be done in the next section. The second term of (^8|a,b) is a single integral, it comes 
about when the boxes act on one internal and on one external diffuson. This corresponds to the case where one 
internal diffuson is almost empty, while the other diffuson contains a arbitrary number of scatterers. In two and three 
dimensions the momentum integral diverges, since for large Q it behaves as J dz J d'^^^QCmtiQ] ^i^) ^ J d'^^^Q Q^^ . 
But when summing the a, b and c contribution this term cancels. The sum gives 



The last two terms involve a double integral describing interference vertices at different points in space; technically 
it arises when both boxes act on external diffusons, or from terms where they do so after partial integrations. This 
term is absent in the expression F^, which contains only one z— dependence as can be seen from Fig. ||(c). When 
performing the integral of F over the transverse momentum Q, the double integral term behaves at large Q as 
J dz J dz' J d'^~^Q£^^^{Q; z, z') ~ J d'^~^Q Q^^ . It is thus convergent. One expects that finally only this term will 
survive. It is also the only contribution depending solely on derivatives of the external diffusons. 

In the quasi one-dimensional case where absorption and frequency terms are absent, the expressions reduce to 



The second part is a well known result for the UCF in one dimension, but a singular part is annoyingly present. 
Before we can obtain the UCF and correlation functions, we have to show its cancellation. We expect the unphysical 
divergency to disappear by summing all diagrams and in the next section we set out finding all leading diagrams. 




(49) 



F,(0) + Ffc(O) + FM = -^mL + 4- 



(50) 
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V. GENERATING DIAGRAMS THROUGH A GENERALIZED KADANOFF-BAYM TECHNIQUE 



Finding the correct and complete set of diagrams by educated guess proves very difficult, especially for the low order 
diagrams. In this section we pursue the technique founded by Kadanoff and Baym | |3^ , |35t to generate all diagrams. 
This method provides a systematic way to construct the diagrams in a particular approximation. The approximation 
is made on the level of a generating functional. The form of this functional is guided by intuition or by prior knowledge 
of the self-energy or some other physical quantity. The theory of Kadanoff and Baym, of which the basics will be 
reviewed in the course of this section, prescribes which diagrams are to be included in the perturbation theory on 
any level in the hierarchy of the many particle Green functions. It was proven by Kadanof and Baym ]34[|3^ that 
this procedure provides a conserving theory in the sense that sum-rules based on conservation of particle number, 
momentum, angular momentum or energy are fulfilled. However, it does not guarantee that a sensible theory in a 
physical or even in mathematical sense will be produced. (The Green functions themselves may not posses the correct 
analytical structure, i.e. they may be non-Herglotz |^^). 

The theory of Kadanoff and Baym is defined for interacting electron systems. It is our goal to describe the scattering 
of light in a disordered medium. The connection between electronic and optical disordered media can be put on an 
solid footing as the equations which describe these phenomena can be mapped onto each other The equivalence 
of interacting and disordered systems, however, is less obvious. Although one can rewrite models of disorder in such a 
way that an effective interaction is present ||37|| , the precise form of that interaction depends on the physical quantities 
that are studied. Furthermore, if one wants to represent the disorder as a two-body, central potential, one is limited to 
special disorder distributions (e.g. Gaussian) or approximations (e.g. second order Born). However, these drawbacks 
do not hinder us at using this theory in a heuristic sense. 

The input of the theory is an approximate functional $ from which we may derive the perturbation theory. In the 
choice we are lead by the theory of weak localization and universal conductance fluctuations (UCF) . The maximally 
crossed diagrams, which are the most important quantum corrections, are generated by the functional depicted in Fig. 
^(a). These wheel-like diagrams are generalized Fock diagrams in the sense that the first term in the sum produces the 
Fock diagram for the self-energy. The UCF diagrams are generated by a functional which may be called " generalized 
Hartree" in the same sense, see Fig. ||(b) and||(c). It turns out that the weak-localization effect and the UCF are 
connected and may be generated from a single diagram |^^. This observation leads to the conjecture that the most 
important contributions in the intensity fluctuations stem from the same diagrams. We thus adopt the functionals as 
given m Fig. |(b) and|(c), as the ones we start our calculation with. 

Let us shortly review the aspects of the Kadanoff and Baym theory we need. Starting with some functional <I> as 
described above, the self-energy SI and the two-particle irreducible vertex El? may be found from 

= ,5$ / SGl (51a) 
= S^^ I (6G\ SGI), (51b) 

where the numbers 1, 1, 2, • • • denote the collection of variables (fci, wi), {k[,uj[), {k2,uj2) ■ ■ ■ and Gl = — (T^'(1)'J'+(1)). 
At this point we generalize these results to obtain irreducible three- and four-particle vertices 

= / (SGl SGI SGl) (52a) 
•Ellll = / (SGl SGI SGl SGt). (52b) 

Note that these generalizations, as the ones that will follow, are not free from internal symmetry. Since we are 
primarily interested in the structure of the diagrams, rather then their multiplicativity, we will make no effort to 
remove this symmetry. The one-particle Green function GJ is connected to the self-energy by the Dyson equation. A 
two-particle Green function GH is connected to the irreducible two-particle vertex by the Bethe-Salpheter equation 

Gl = Gl{0) + Gl (0) Gl (53a) 

12 — <-ri <-r2 + "-^3 ^34 "-^42 1 [OM) 

where GJ(0) indicates the bare propagator. Integration over the free variables occurring as super and subscript is 
assumed. Again, we may generalize these results to obtain a three-particle Green function GUI ^^"^ ^ four-particle 
Green function GJgJ. We find || 

/^I23 //^I /^53 I /^I3\ / rS r2 I V45 /^24 \ , /^I /^4 / ^45 /^235 i V'456 /^38 /^25N /^An\ 

(-^123 — V^5 (-^13 + '-"l '~'53)["2 "s + ^45 <-^24 ) + ^i^l [^li ^235 + ^456 '"^36 '"^25^ (,04aj 

GI234 //^I34 /^B I /^I3 /^54 , /^I4 /^53 i /^3 /^534'\ /• r5 r2 , ^56/^26^ i 

1234 — l'-^B34 '-^1 + '-"53 '-"l4 + '"^54 '-"13 + 5 '-^134Jl"2 + ^56 '^26-' + 

(G;« Gl + Gl GJ«)(E5?I GZ Gil + nt GllDiSl 61 61 51 + SI S'l <5| S'l) + 

1^1 /v^STS /^4S /^37 ^26 i ^HBl ( /^3i7 /^26 i /^37 ^246 i /^47/^236'. i ^56 /^2346\ ('^A\^\ 

'^b "-"l l-^BSTS '-"48 '-"37 '-"26 + ^567 l'-"347 '""26 '"^37 '-"246 '-"47'-"236J ^56 '-"2346 j 5 {OW ) 
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where the 5's are Kronecker delta's. These equations provide the full information on the fluctuations in transport 
quantities given the exact set of irreducible vertices , Sj;? , Sj?? , An approximation of this set can be found 

by choosing a suitable functional $. 
The functional we adopt is 



$ := Gl Gl -SI 



E 



(55a) 



where ii5[Ti] and Cj|[?^] are the product of n scattering terms in the particle- hole and particle-particle channel 
respectively. The last two terms are the ones depicted in Fig.^(b) and ||(c). Algebraicly we have 



5^1 = Li? [l]=Cll[l]^5lT5l 
L\l [n] = S'jI Gl Gl 51: • • • S-^ 

C12 r ] C34 C 044 

Tnl'll — Otk (_r. (j-^Ojs • • ■ J-- , 



(56a) 
(56b) 
(56c) 



where T is the scattering strength in second order Born approximation. After differentiation of $ with respect to G 
ladder diagrams and maximally crossed diagrams appear. Each element in these series have prefactor one. Therefore 
we define the full ladders in both channels as 



-^12 — ^11 + ^13 G\ Gl Lll — ^ Ljj [n] 



G3 /^12 
4 *-^3 ^43 



OO 

E 

n=l 



Gl 



(57a) 
(57b) 



Systematic application of Eqs.(pl[) and (52) on the functional ( |55| ) leads to the following irreducible n-particle 
vertices 



■srii2 

^12 
V-*123 

^123 



y>1234 
^1234 



2 1^21 ^21 ~ ^2l) 

L12 I /^12 Q12 I /^3 r^4 / T 32 T 14 , T 34 T 12 i r^32 /^14 i r^23 /^14 \ 

21 "T *-^21 ~ '-'21 "T '-^3 ^4 \^14 ^32 "T ^12 ^34 Ogg "h Oj^ Ogg J 

PerS {G5(L-L- + C?|C-) + 

r^i r^l 7- 16 7- 57 7- 72 I 7- 16 7- 57 7- 32 I , f~ll6 /^73 /^62 , ,^16 /^37 ,^52 N \ 

L75 LXg (J-^ I, 3^52 -'^13 + ^52 -'^13 -^^76 3 "-^52 "-"IB '""37 + "-^52 "-^16 '-'73 ) J 



G5 /^6 / r 35 r 45 r 12 
5 '-'6 1-^14 -'^62 

G5 /^6 /^7 /^S /" , 7-54 r S3 7- 72 r 
5 '-'6 '-^7 '-^sU-^lS -^47-^36 ^ 
1 /^84 >^73 /^52 >^16 , x^48 ^73 y^52 y^l 

s'-^Tfi "-^is "-^s? "-^sq "t '-'tb '-'aa "-^57 "-^s^ 



, 7-34 7-65 7-12 I , 7- 35 7- 16 7- 42 
5-'^I8 -'^42 -'^35 + 7-^^12 -'^34 -'^BS ^ 



^35 

LI 



I /^46 /^35 /^12 
r "-^H '-'B4 '-^SB " 



/^64 /^35 /^12\ I 
'-'IS "-"iB '-"35-' "T 



1 >^48 >^37 >^52 y-tl 
'7*— 'Tr '—'5^ 



1 /^37 <^52 <^84 /^16\\ 
I'-^iS '-'73 '-'52 '-'48 Jj J 



(58a) 
(58b) 

(58c) 



(58d) 



where the operator Perj?; ;^ produces all the permutations of its operand in the variables (1, 1), (2, 2), • • • , {n, n), i.e. 
nl diagrams. Expansion of Eq.(|5^) using Eq. (|58|) is obviously a considerable task. Taking into account that we only 
need those diagrams which are different in a topological sense, the effort stays manageable. As usual, diagrams which 
contain loops do not contribute. 

To illustrate the procedure set out above we have a closer look at two examples. First of all, part of the integral 
equation for GHH is needed. The relevant part is depicted in Fig. |^. Then, in Fig. |5|(a) we start off with the functional 
as depicted in Fig. |^(b). Successive functional differentiation (cutting) as prescribed in Eqs.([5l|) and ( |52| ) produce the 
one-, two-, three- and four-particle irreducible vertices, given in Figs||(b), ^(c), ||(d) and ||(e) respectively. Insertion 
(glueing) of the latter irreducible four-particle vertex in Fig. ^ amounts to adding external diffusons. Of course, more 
diagrams are produced in this particular example but for the sake of simplicity no attention is paid to these diagrams. 

In the second example we start off with the diagram containing the maximally crossed vertex. Figs, ^(a) or 0(c). 
By functional differentiation, see Figs ||(b) and ^(c) and subsequent integration in the Bethe-Salpheter equation (|53|) 
we obtain a contribution which consists of two maximally crossed diagrams with a dressing. This element is the input 
for Eq.(54b), see Fig. ^, and generates an important contribution to the whole set. 

The set of diagrams thus generated by the Kadanoff-Baym approach was verified using a computer program. We 
developed this program to generate all diagrams with six or less scatterers automatically. After generating the set, 
the program checks it against double counting. Next, momenta are assigned to all propagators such that momentum 
conservation is obeyed. It then determines whether a particular diagram is of leading order in (k£)~^; this is the case 
if all propagators can have a momentum approximately equal to k. The leading diagrams are expressed in terms of 
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the standard integrals defined in section Vl. The subsequent analytical calculation and summation of diagrams was 
done by hand. The program turns out to be especially useful in determining the precise set of leading short-distance 
diagrams, but also the long distance diagrams of section III were reproduced. It should b e men tioned that apart from 
the diagrams from the Kadanoff-Baym approach we also found extra diagrams; in section VI C they are discussed and 
it is shown that they cancel. 



VI. CANCELLATION OF SHORT DISTANCE CONTRIBUTIONS 



In this section we show that the strong divergence, that occurs if both Hikami boxes act on internal diffusons, 
cancels. Since a diffuson decays rapidly if its momentum becomes large, terms of the diffusons with few scatterers are 
dominant when the momentum is large, which cause the divergence. The cancellation is thus shown by considering 
in great detail the short distance processes. 

It turned out that complications arise for diagrams with less than four scatterers. The external diffusons may still 
contain an arbitrary number of scatterers. For simplicity the calculation in this section is performed in the bulk and 
the external diffusons are not attached to the diagrams. Since the divergency is independent of the external momenta, 
the cancellation is generally proven at zero external momentum. For simplicity we may then consider an infinite 
system. As a result all factors can be expressed in the d-dimensional momentum q = {Q,qz)- We first look at fixed 
value of the d-dimensional internal momentum q and postpone the integration. We will show that under present 
conditions the integrand is zero for all q, so that there is no divergency after integration. 

All diagrams can be factorized in products of the integrals 

4"r = ^ / (^Gnp)G*'(p)G™(p + q)G™(p + q). (59) 

In the calculation of the diffuson and the Hikami vertices in the previous sections we expanded the integrals in q£, 
Since we are after contributions for q ^ 1/^, this expansion is not allowed. We can still assume that q <C fc, as we do 
not need the physics on length scales comparable to the wavelength, but comparable to one mean free path only. The 
integrals needed in the calculation are given in Table ^, where Ai is an angular average defined as 

A.(,) = ^/#(^^^ (. = 1,2,3). (60) 

The integral for i = 1 yields the diffuson kernel Ai = 'Ai-ct'An{q£)/q£. The internal diffuson now reads 

A„.(.) = f^. (61) 

We recover previous results by expansion in {q£). This yield Cint = 127r/(£'^q^), agrees with the bulk solution of (p2). 



A. Short distance properties of the long range diagrams 



First, look at the long range diagrams with few scatterings. By writing out the dressing of the vertices one obtains 
the detailed structure of the diagrams, given in Fig. ^ In the diagrams Fig. |(a)i-(c)i, the incoming external diffusons 
are connected to 1 and 3. The diagrams in Fig. |^(a)ii-(c)ii are obtained by connecting external diffusons 1 and 2 to the 
incoming side. If such a diagram contains only a few scatterers, it can be grouped into one of the classes of the long 
range diagrams. A simple example is the upper left diagram of Fig. 0, drawn there with six scatterers. Interchanging 
the external diffusons numbered 2 and 3 clearly leads to another topology. However, with only two scatterers present, 
the topology does not change under this operation, and one must be careful not to overcount this term. Bearing this 
important observation in mind, we sum all long range diagrams for an arbitrary number of scatterers. The expression 
for each diagram and its combinatorical factor is given in Table ^ It is convenient to collect diagrams with an equal 
number m of scatterers that connect different propagators. Scatterers on which one given amplitude scatters twice 
are thus, momentarily, not counted as new scatterers. Looking at Fig. |^, m equals the total number of scatterers 
minus the number of scatterers indicated with curved lines. After this resummation we find the important result that 
each class of "long range" terms vanishes for m not equal to 2. Denoting the terms of Fig. respectively by Ra, Rb, 
and Rc we find Ra = Rb = Rc = '^i m > 2. If m = 2, the classes cancel against each other, since we then obtain 
Ra = Rb = -Rc/2 = -£tAI. This imphes 
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Ra+Rb + Rc = 0. (62) 

We have thus shown that the long range diagrams for an infinite system cancel and thus cause no divergency after 
integration over q. It is essential that the degeneracy of the low order contributions be counted properly. 

Another, more standard way to verify the important cancellation is the following. Recall that we look at a bulk 
situation with fixed internal momentum q, while the external momenta are put equal zero. Beyond the diffusion 
approximation an internal diffuson is given by Eq.(|6T|). We also need the Hikami-boxes beyond the diffusion approx- 
imation. It holds that 

i/4(q,0,-q,0) = llUq) + ^ {/i?(<z)}' + y {I'^Hq)}' ■ (63) 
Inserting the values from the table and doing similar but longer calculations for two other vertices we find 



i74(0,q,0,-q) = 2i/4(q,-q,0,0) - ^Ai(l- Ai), (64a) 



^6(0, 0, q, 0, 0, -q) = Y^^i(l - ^i)(l - 3Ai). (64b) 
Now the long range diagrams of Fig. ^ can be evaluated. For the internal part one has at fixed q 

2 Hi{q, 0, -q, OfClM + iH^q, -q, 0, 0)^ClJq) + ^H^iq, 0, 0, -q, 0, 0)A„t(g) 

= ^{2A? +AI + 2Ai(l - 3Ai)} = ^(2yli - iAl). (65) 

It is essential that the denominators 1 — Ai have disappeared from this expression. As mentioned, it means that all 
high order terms cancel, allowing a cancellation of the remainder by low order contributions. From Tabic || it is seen 
that we have overcounted six types of low order terms. The correction to be subtracted is 

— { 4 jiAi + ^2 + ^3) 



^1 +8:^:^(2^1+^2)^^1+8— —-(A2 + 2A3) 



)V87rfc2 327rfc3' ' 87rfc STrfc 327rfc3 

It indeed exactly cancels expression (|6^), which arose as a remainder of the long range terms. 

From the results of present section one may be tempted to conclude that the long range diagrams do not lead to 
divergencies if the correct degeneracies of the low order terms is properly taken into account. Though this conclusion 
is correct, it is too early to draw it, since we show in next parts that some further complications arise. 



B. Extra short distance contributions 



Unfortunately, we have not yet finished with the calculation as also other terms are of leading order at short 
distances. First, diagrams with a different topology also occur. There are five new diagrams, all with two scatterers, 
see Fig. ^. They can be looked upon as diagrams from the set of Fig. 0(a), but without any scatterers in the ladders. 

Second, diagrams can have more than one configuration for a resonant arrangement of the momenta, i.e. several 
configurations can be leading. The most important contribution to the integral (59) arises if the amplitude loop 
momentum p is close to the wavevector |p| = k. With more loops present it is sometimes possible to find more than 
one choice for the momenta of the amplitude propagators. This is illustrated in Fig. ^ where the diagram of Fig. 0(a)l 
is repeated. The expression for the shown diagram has contribution from three different poles. The first arrangement 
of the momenta is dominant for any number of scatterers, this one is part of the set of long range diagrams for 
any number of scatterers. The two configurations on the right give an extra leading contribution. This only occurs 
with two scatterers present; with more scatterers the two diagrams on the right become subleading, and only the left 
diagram remains. Apart from this diagram, the same thing occurs for the diagrams of Fig. |^(c)l and Fig. 0(c)2 also 
with two scatterers, and for the diagram Fig. |^(c)6 with three scatterers. 
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Also the diagrams of Fig. || have more than one pole: The four left diagrams of Fig. ^ have two resonant momenta 
configurations, the right one has three. We sum all the contributions thus found (41 in total) and denote them as Sn 
for n scatterers. They are only non-zero for diagrams with two or three scatterers, and read 



^2 
^3 



An 



jQ.l j0,2 j^O.O 
^1,1 ''1,1 ^1,2 



2fc4 



A1A2 



One sees that these special contributions cancel. 



2k 



-AIA2 



(67a) 
(67b) 



C. Interference vertices without partner exchange 



Apart from the diagrams generated by the Kadanoff-Baym approach, we found another class of leading long range 
diagrams when we generated the diagrams by computer. In these diagrams the two diffusons have two scatterers in 
common, but no amplitude is exchanged. They either have one internal diffuson, see Fig. |l^ or two, see Fig. These 
diagrams are actually of similar type as Fig. g. The diagrams with one internal diffuson is analogous to Fig. ||(c)i, 
while the ones with two internal diffusons are analogous to Figs. ^a)i,(b)i. The equivalents of Figs. ^a)ii,(b)ii and 
(c)ii also occur. The important difference with the previous diagrams, though, is that at the interference vertices 
no amplitudes are exchanged. In this respect they differ qualitatively from Hikami-boxes, where partner exchange 
does occur. Explicit calculation shows that these terms cancel both at zero and non-zero external momentum, this 
was also already noted by Kane, Serota and Lee |23|, see their Fig. 5(c). These classes of diagrams can thus be fully 
neglected. This cancellation is due to time reversal invariance. Apart from these two classes, the Kadanoff-Baym 
approach, using (^5|), generates all presented diagrams as we checked with the computer program. The generating 
functional for the extra diagrams of Fig. |l^ is drawn in Fig. |ll|. 

In conclusion, the analysis of this section VI shows that when the external momenta vanish, all leading terms 
cancel at fixed value of the internal momentum. Upon integration over the loop momentum one still has zero and, 
in particular, not a divergent contribution. For an infinite system the theory is thus well-behaved at short distances. 
The divergence has canceled in a careful study of the short-distance process; renormalization as known from field 
theories was not needed. 



VII. APPLICATION OF CONDUCTANCE FLUCTUATIONS TO OPTICAL SYSTEMS 



We have now seen that in infinite systems all divergent terms cancel. In realistic systems, such as a slab, the same 
diagrams describe the relevant physics. They have to be evaluated with appropriate diffuson propagators. Generally 
they can be written as the bulk term with additional mirror terms. For a slab they were given in section ||. Knowing 
that the short distance behavior is regular, we can selfconsistently consider all scattering diagrams in the diffusion 
approximation. For the long range diagrams this was done already in section IV, where all long range contributions 
were evaluated. The discussion of previous section has shown that the only new effect comes from the subtraction of 
Eq. (|66|). In the diffusion limit for a quasi one-dimensional system at fixed transversal momentum Q it results in a 
contact term, labeled Fd- 



Fd{M) = / dzCU^) 



2tt 4fc4 



4fc4 



m 



(z). 



(68) 



This indeed cancels exactly the leading divergency that remained in Eq.(^) for the sum of Fa, F^ and Fc. The milder 
divergency present in the individual terms Fa, Fb, Fc in two and three dimensions canceled already by summing 
them. Diagrammatically this can be seen as follows: The diagrams responsible for the latter divergence contain one 
internal diffuson (such in Fig. ^(c)). For these diagrams one does not have the complications of the previous section; 
double counting corrections, different resonant momentum configurations, and extra diagrams are absent. This was 
confirmed by our computer generated diagrams. The only contributions which remain in ( p9[ ) are thus the double 
z— integrals. This finite remainder of our theory yields the sought conductance fluctuations. We thus obtain, with 
again /i4 = f.^ / (ASTrk'^) 
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F{Q, K, n) = Fa{M) + F,{M) + F,{M) + Fd{M) 

= Ahl [ /dzdz'A„t(z,z';Af)A„i(z,z';M*) [/:£(^)£^2^^(/) + 



+ hlj J dzdz' [CUz,z';M)+Cl,{z,z';M*)] ^ [CUz)Cout{z)] ^ [£in(z')^out(^')] , (69) 

in which again Af^ = + + iil. This equation is the central resuh of the present paper. The upper hne of ( |69| ) 
corresponds to the diagrams of Fig. |^(a), whereas the lower line corresponds to the diagrams of Fig. |(b). Note that 
only derivatives of external diffusons are present. As compared to the first line, extra terms are present in the second 
one. According to (^2|) we have 

^ [CUz)C,utiz)] - 2£;„(z)£;,,(z) + 2K'C,^{z)C,^,iz) (70) 



The K— terms are extra terms arising when absorption is present. Finally, with Eq.(38c) the value at vanishing 
transversal momentum gives the variance of the conductance in one dimensions, integration over the transversal 
momentum yields the correlation in tw o an d three dimensions. 



Using the general result of ( |69| ) and ( 38c ) various cases are considered by inserting the diffusons derived in section 
||. First consider the case of fully transmitting surfaces; if we neglect absorption and frequency differences this gives 



3 2 + 2Q^L^ - 2 cosh 2QL + QL sinh 2QL 

2 Q^U^^^FoT 



which decays for large Q a,s Q In this case we recover 

(T^)^ ^ — fn 0.133, quasi Id (72) 
15 

3 W W 
= ^C(3)— ~ 0.116—, quasi 2d (73) 

TT-^ L L 

1 ^2 

-^7^-0.159^ 3d, (74) 

in which ^ is Riemann's zeta function. These are well known results Q|. We determine also the frequency dependency 
of the correlation; this is important as it determines the frequency range of the light needed to see the fluctuations. 
Taking the frequency dependency into account we obtain 

4 (M*2 - ]VP + M^M*LcothM*L - M*^MLcothML) 
FiQ,Lo)^^ ' 



*2^ 



2 + 2APL^ - 2 cosh 2ML + ML sinh 2ML , , 

+Rc — 5 . (75) 

2M^L^ sinh^ ML 

The correlation decays for large frequency differences as fi^"''/^, as was stated by Lee, Stone and Fukuyama 

We were unable to perform the double integral over the position analytically in presence of absorption. In Figs. |l^, 
|l|, and |l| we show the one, two, and three dimensional correlation functions for various values of the absorption. It 
is seen that especially the top of the correlation is reduced due to absorption. 

Next we applied our theory to the case of partial reflection at the surfaces of the sample. We assume an index 
of refraction m — \/ ep / ei ^ 1. For our purpose the internal reflections are coded in only one parameter, the 
injection depth zp, see (|l8|), (p6|), (|35|). For optically thick samples the correlations are now determined by zp/L. It 
is physically clear that the internal reflections lead to a less steep diffuse intensity in the sample as function of the 
depth. The fluctuations are proportional to the space derivatives and are thus reduced. The results are presented in 
the Figs. |l^, and 18 where the correlation functions is plotted for various values of the ratio between extrapolation 



length and sample thickness. One sees that already for small ratios of zp/L the correlation is significantly lower than 
without surface internal reflection. It is seen that neither the variance (the value at vanishing 17), nor the form of the 
correlations are universal, they are sensitive to absorption and internal reflections. 
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VIII. DISCUSSION AND OUTLOOK 



The universal conductance fluctuations (UCF) of mesoscopic electronic systems have a direct counterpart in other 
mesoscopic systems with multiply scattered classical waves. The corresponding average normalized correlation function 
is the C3. It is known that a naive calculation of this object in the Landauer approach is plagued with short 
distance divergencies. In this work we have presented a detailed diagrammatic approach to the calculation of the C3 
correlation function. We first evaluated the leading long range diagrams. As expected, each diagram, but also their 
sum, contains a short distance divergency. Also a subleading divergency occurs, but for a slab geometry is was found 
to cancel automatically. The study of the cancellation of the leading divergency then was the main theme of the 
paper. Consistency of the approach requires finding extra contributions that exactly cancel the already determined 
divergency, and therefore serves as proof that the set of long range diagrams exhaust all of them. 

We developed a diagrammatic method that systematically generates all leading scattering diagrams. This set was 
checked by computer. A large number of diagrams had to be considered in detail. We summed these diagrams for an 
infinite system at fixed value of the loop momentum. It was seen that also beyond the diffusion approximation the 
action of Hikami boxes on the long range propagators is to eliminate the long range terms, and leave only some low 
order contributions. Moreover, some low order terms of diagrams with diffusons have a lower degeneracy than their 
higher order equivalents. Taking this into account led exactly to a cancellation of all terms in an infinite system. Next 
we discussed that some extra classes of leading diagrams occur, but they all add up to zero. Thus in an infinite system 
all diagrams cancel, so that, in particular, no short range divergency occurs. All short range contributions could be 
coded in a contact term (|68|). It would be interesting to investigate how it is derived in a non- linear sigma-model 
formulation of the theory, possibly along the lines of Serota, Desposito and Ma . 

Subsequently the theory was applied to systems of finite size. Here the short distance divergencies cancel as well, 
because the large scale geometry of the system does not have infiuence on short range effects. The final result is non- 
zero as it describes the correlation function in terms of derivatives of external diffusion propagators of the geometry 
considered; such terms have no meaning in an infinite system. 

Our central result for the correlation function of the conductance is given in (|69|). It is obtained by adding the result 
of the long range diagrams (|4^) and the contact term (^8|). When there is no absorption, agreement is found with 
the result of Kane, Serota and Lee All external diffusons are differentiated once. When absorption is present, 
however, their approach is no longer valid. We found that extra terms appear where some external diffusons are 
differentiated twice or are proportional to k^, see (|70|). 

We have applied the results to realistic optical systems. The frequency dependent C3 correlation function was 
calculated for the case where a diffuse incoming beam is used and all outgoing intensity is collected. It was seen that 
both absorption and internal refiections decrease the correlations by a considerable amount. This is important for a 
quantitative analysis of experimental data. 

Electromagnetic measurements that involve the C3 correlation have been reported by Genack et al. |^]. These 
authors were able to describe the data of their infrared experiments by adding the Ci, C2 and C3 contributions, but 
they incorrectly assumed that the C3 is frequency independent. The experimental investigation of optical universal 
conductance fluctuations is known to be very difficult. One problem is that if the incoming beam has to be diffusive, 
it will have a low intensity. 

We propose here a different way to measure the same interference effect. Consider a laser beam coming in in a 
given direction a and measure the frequency dependent total transmission. Such can be done using an integrating 
sphere . Then repeat the measurement for a very different incoming direction c. Each of these two signals will 
exhibit the large C2 correlation function However, when the directions a and c are not close to each other, the 
C2 will not contribute to the cross correlate. The cross correlation of the total transmission is much smaller than the 
autocorrelates; it just represents the typical term of the UCF in (|39|). As this cross correlate is of relative order 1/5^, 
it is of the same order of magnitude as the third cumulant of the total transmission. A very precise measurement of 
that quantity was carried out by de Boer, van Albada and Lagendijk, and reported in collaboration with two of us 
[0. It may therefore be expected that it is just possible to measure C3, and thus essentially the UCF, with visible 
light. 
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file=c3figl .eps,width=8cm 

FIG. 1. Left: an example of an actual scattering process; a retarded (full line) and an advanced amplitude (dashed line) 
come from the left and share the same path through the sample. Right: schematic representation of the average process, the 
diffuson. 



file=c3fig2.eps,width=16cm 

FIG. 2. The leading contributions to the conductance fluctuations, apart from some special short distance processes dealt 
with in section The incoming diffusons from the left interfere twice before they go out on the right. The close parallel lines 
correspond to diffusons; the shaded boxes are Hikami vertices; q denotes the free momentum which is to be integrated over. 

file=c3fig3.eps,width=8cm 

FIG. 3. The generating functional which produces diagrams for (a) weak localization, (b) and (c) the UCF. 
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FIG. 4. Part of the integral equation from Eq.(54b). First indicated term on the r.h.s. produces the external ladders. 
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FIG. 5. Particular example how the double diffuson structure is produced. Functional differention with respect to G is 
indicated by cutting out a line (pair of scissors) from (a) to (e). Pasting the result in the integral equation, indicated in Fig. 
^ results in adding external ladders (glue) . 
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FIG. 6. Particular example how a set of cooperons is produced. Functional differention with respect to G is indicated by 
cutting out a line (pair of scissors) from (a) to (e). Pasting the result in the integral equation, indicated in Fig. ^ produces a 
four particle quantitiy (and adds external ladders). 
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FIG. 7. The dressings of the shaded vertices of Fig. H are written out explicitely, this shows the detailed structure of the 
long range contributions to the conductance fluctuations. The close parallel lines, only explicitly shown in the first diagram, 
are the external diffusons. Fig. |^a)i corresponds to the case were 1 and 3 are incoming diffusons; Fig. H(a)ii corresponds to 
the case were I and 2 are incoming diffusons. The vertical, diagonal or curved lines linking two dots represent a common 
scatterering of two amplitudes. From top to bottom the horizontal lines in a diagram are advanced, retarded, advanced and, 
retarded propagators, respectively. The number of scatterers in the internal and external diffusons is arbitrary. Topological 
equivalent diagrams are not shown. 
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FIG. 8. Leading diagrams with topologies not yet contained in the long range diagrams. They all contain two scatterers. 
file=c3fig9.eps,width=8cm 

FIG. 9. Different choices for the loop momenta, all yielding leading contributions; thick and thin lines depict different large 
momenta, both with length k. The first diagram is leading for any number of scatterers. When only two scatterers occur also 
the other two diagrams are leading. 
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FIG. 10. Leading diagrams without amplitude exchange with one internal diffuson. The two drawn diagrams cancel against 
each other. All diagrams of this class cancel. 
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FIG. 11. The generating functional that has to be included in the Kadanof-Baym approach to generate the diagrams of 
Fig. 0. 
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FIG. 12. Leading diagrams without amplitude exchange with two internal diffusons. The two drawn diagrams cancel against 
each other. All diagrams of this class cancel. 
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FIG. 13. The correlation of the conductance as a function of the frequency for various absorption strenghts in one dimension 
without internal reflection. From upper to lower curve: no absorption(«; = 0), k = l/L and, k = 2/L. 
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FIG. 14. The correlation of the conductance as a function of the frequency for various absorption strenghts in a two dimen- 
sional slab without internal reflection. From upper to lower curve: k = 0, k = 1/L and, k = 2/L. 
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FIG. 15. The conductance correlation function versus the frequency difference for various absorption strenghts in a three 
dimensional slab without internal reflection. From upper to lower curve: k = 0, k = 1/L and, k = 2/L. 
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FIG. 16. The influence of internal reflection on the one dimensional frequency correllation. With zo = 0,L/10, L/5 (upper, 
middle, lower curve); no absorption. Also here the fluctuations are reduced. 
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FIG. 17. Influence of internal reflection on the frequency correllation in 2D; zq — 0, L/10, L/5 (upper, middle, lower curve); 
no absorption. 
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FIG. 18. Influence of internal reflection on the frequency correllation in 3D; zo — 0, L/10, L/5 (upper, middle, lower curve); 
no absorption. 
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TABLE I. The short distance diagrams can factorized in these integrals, defined in Eq.(|59|); Aj is defined in the text. 
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TABLE If. Table used in the calculation of the diagrams presented on Fig. ^ for zero external momenta. The expression 
for each diagram is given in the second column. Its combinatorical prefactor in the other columns for different number of 
scatterers. The last column counts the number of scatterers that are not included as such in the resummmation. There are six 
diagrams with degeneracy lo-wer than expected: cl for n — 1, al,c2,c5 for n = 2, and c6,c7 for n = 3. 
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